RIFFLE SHUFFLES OF A DECK WITH REPEATED CARDS 



SAMI ASSAF, PERSI DIACONIS, AND K. SOUNDARARAJAN 

Abstract. We study the Gilbert-Shannon- Reeds model for riffle shuffles and ask 'How 
many times must a deck of cards be shuffled for the deck to be in close to random order?'. 
In 1992, Bayer and Diaconis gave a solution which gives exact and asymptotic results 
for all decks of practical interest, e.g. a deck of 52 cards. But what if one only cares 
about the colors of the cards or disregards the suits focusing solely on the ranks? More 
generally, how does the rate of convergence of a Markov chain change if we are interested 
in only certain features? Our exploration of this problem takes us through random walks 
on groups and their cosets, discovering along the way exact formulas leading to interesting 
combinatorics, an 'amazing matrix', and new analytic methods which produce a completely 
general asymptotic solution that is remarkable accurate. 



1. Introduction 

A basic question in scientific computing is 'How many times must an iterative procedure 
be run?'. A basic answer is 'It depends.'. In this paper we study the mixing properties of 
the Gilbert-Shannon- Reeds model IT3] for rime shuffling a deck of n cards and ask how 
many times the deck must be shuffled for the cards to be in close to random order. Our 
answer depends not only on the metric we use to measure distance to uniformity, but also 
on the particular properties of the deck that are of interest. 

To be precise, we consider a 'deck' to be a multiset of n cards. We shuffle the deck by 
first cutting it into two piles according to the binomial distribution, and then riming the 
piles together by successively dropping cards from either pile with probability proportional 
to the size. This process defines a measure, denoted Qi{a), on the symmetric group S n . 
Repeated shuffles are defined by convolution powers 

(1) Qf(a)= Y, Q2(r)Qf- 1) (u J ). 

LU-T=(T 

This shuffling model, which accurately models how most people actually shuffle a deck 
of cards, was introduced by Gilbert and Shannon [Tl] and independently by Reeds [13] . 

Bayer and Diaconis [3] generalized this to o-shuffles, which is the natural extension to 
shuffling with a hands: the deck is cut into a packets by multinomial distribution and cards 
are successively dropped from packets with probability proportional to packet size. Letting 
Q a (o~) denote this measure, they show that convolution of general a-shuffles is as nice as 
possible, namely 

(2) Qa*Qb = Qab- 

Thus it is enough to study a single a-shuffle of the deck. 

To that end, denote the uniform distribution by U = U(o~). For a deck with n distinct 
cards, U = 1/n!, and for a more general deck with D\ l's, D% 2's, up to D m m's, we have 
^ = V(^i a!" 1 )' There are several ways to measure the distance between Q a and U, 
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though for the purposes of this paper we restrict our attention to total variation distance 
and separation distance. 

The total variation distance is defined by 

(3) \\Q a -U\\ TV = max \Q a {A) - U{A)\ = \Y j \Q a {a) - U{a)\. 

subsetsA I z — ' 

a 

In general, the formulas for Q a (o~) may be quite complicated, making calculations of total 
variation intractable. Therefore we will also consider the separation distance defined by 

(4) S EP(a)=maxl-^4. 

c U (crj 

Here, only a single probability needs to be computed, though as we shall see even that 
can be difficult. From the formulas above, one can easily see that separation provides an 
upper bound for total variation, which makes separation a good measure to use when total 
variation becomes too complicated to compute. 

In widely cited works, Aldous [1] and Bayer and Diaconis [3] show that |log 2 (n) + 
c shuffles are necessary and sufficient to make the total variation distance small, while 
21og 2 (ra) + c shuffles are necessary and sufficient to make separation small. These results, 
however, look at all aspects of a permutation, i.e. consider a deck with distinct cards. 
In many card games, only certain aspects of the permutation matter. For instance, in 
Baccarat, suits are irrelevant and all 10's and picture cards are equivalent, and in ESP card 
guessing experiments, a Zener deck of 25 cards with each of 5 symbols repeated five times 
is used. It is natural, therefore, to ask how many shuffles are required in these situations, 
and so we consider a deck to have repeated cards. 

Many results are known for how long it takes certain features of a permutation, e.g. 
longest cycle, descent structure, etc, to become random; for a thorough treatment of such 
results see [8]. The particular problem we address in this paper was first addressed by 
Conger and Viswanath (5J [7] who derive remarkable numerical procedures giving useful 
answers for cases of practical interest. 

In this paper, we present many of our main results from [2], giving exact formulae and 
asymptotics for a deck of n cards with D\ cards labelled 1, D2 cards labelled 2, . . ., D rn 
cards labelled m. Our results are proved from the deck starting 'in order', i.e. with l's 
on top through m's at the bottom. In Section [2j we show that the processes we study are 
Markov by framing the problem in the context of random walks on cosets. We derive a 
formula for the transition matrix following a single card in Section [31 and show that this 
matrix shares many properties with Holte's 'Amazing Matrix' [12]. In Section^ we consider 
a general deck, limiting our metric to the separation distance, and derive new formulae and 
asymptotic approximations which we unify into our 'rule of thumb' formula. Section [5] 
shows that our results depend on the initial configuration of the deck, a fact also observed 
by Conger and Viswanath [6] . This extended abstract contains precise statements of 
our main results along with the main ideas of the proofs; for full details see [2]. 



2. Random walks on Young subgroups 

In this section, we reformulate shuffling in terms of random walks on a finite group, so 
that our investigation of particular properties of a deck becomes a quotient walk on Young 
subgroups of S n . 
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Let G be a finite group, and let Q be a probability on G, i.e. Q(g) > and ]C g eG Q(d) = 
1. Take a random walk on G by repeatedly choosing elements independently from G with 
probability Q, say 51,52,53, ■ ■ •, and, beginning with the identity element 1q, multiply on 
the left by gi. This generates the following sequence of elements, the left walk, 

1g, 91, 5251, 535251, 

By inspection, the chance that the walk is at g after k steps is given by convolution formula 
© Q* k (g), where Q°(g) = S la>g . 

To focus on certain aspects of the walk, we choose a subgroup and consider the quotient 
walk as follows. Let H < G be a subgroup of G, and let X denote the set of left cosets of 
H in G, i.e. X = G/H = {xH}. The quotient walk on X is derived from the left walk on 
G by reporting the coset to which the current position of the walk belongs. This defines a 
Markov chain on X with transition matrix given by 

(5) K(x, y) = QiyHx- 1 ) = £ Qiyhx^ 1 ). 

heH 

Note that K is well-defined (i.e. independent of the choice of coset representatives) and 
doubly stochastic. Thus the uniform distribution on I, P = is a stationary 

distribution for K. The following result, showing that powers of K correspond precisely to 
convolving and taking cosets, is intuitively obvious with a straightforward proof. 

Proposition 2.1. For Q a probability distribution on a finite group G and K as defined in 
we have 

K l (x,y) = Q* l (yHx- 1 ). 

We may identify permutations in S n with arrangements of a deck of n cards by setting 
a(i) to be the label of the card at position i from the top. For instance, the permutation 
2 1 4 3 is associated with four cards where "2" is on top, followed by "1", followed by "4", 
and finally "3" is on the bottom. Therefore the random walk on S n with probability Qi 
corresponds precisely to riffle shuffles of a deck of n distinct cards. If we consider the first 
D\ cards to be labelled "1", the next D2 cards to be labelled "2", and so on up to the 
last D m cards labelled "m" , then this corresponds precisely to the coset space of a Young 
subgroup, 

X = S n / (S Dl x S D2 x • • • x S Dm ) . 

Thus Proposition 12.11 shows that the processes studied in the body of this paper are Markov 
chains. 

3. A NEW 'amazing' matrix 

Suppose the ace of spades is on the bottom of a deck of n cards. How many shuffles does 
it take until this one card is close to uniformly distributed on {1, 2, . . . , n}? We analyze this 
problem by writing down the transition matrix following a single card through an otherwise 
indistinguishable deck. 

Proposition 3.1. Let P a (i,j) be the chance that the card at position i moves to position j 
after an a-shuffle. For 1 < i, j < n, P a (i,j) is given by 

k=l r=l K J V 7 
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where r ranges from I = max(0, (i + j) — (n + 1)) to u = min(i — l,j — 1). 

Proof. Consider the number of ways that an inverse a-shume can bring the card at position j 
to position i. First, the card at position j must have come from some pile, say k, 1 < k < a. 
Say r of the cards above this came from piles 1 to k, and so the remaining j — 1 — r came 
from piles k + 1 to a. Those r cards all must appear before the card at position j in P'~ ) 
ways. This leaves i — l — r cards below position j which came from piles 1 to k — 1 in f ■ ™~jM 
ways, and the remaining cards must be from piles k to a. □ 

For example, the n x n transition matrices for n = 2, 3 are given below. 

, , /(a + l)(2o + l) 2(a 2 -l) (a-l)(2a-l) 

/ " r ' " ' ^ f 2(a 2 - 1) 2(a 2 + 2) 2(a 2 - 1) 

(a-l)(2a-l) 2(a 2 -l) (a + l)(2a + l) 

These matrices share many properties, given in Proposition l3.2l with the 'amazing matrix' 
discovered by Holte [T2] in his study of the 'carries process' of ordinary addition. Diaconis 
and Fulman [9] show that Holte's matrix is also the transition matrix for the number of 
descents in repeated a-shufhes. We have not been able to find a closer connection between 
the two matrices. 

Proposition 3.2. The transition matrices following a single card have the following prop- 
erties: 

(1) they are cross-symmetric, i.e. P a (i,j) = P a {n — % + 1, n — j + 1); 

(2) they are multiplicative, i.e. P a ■ Pb = P a bi 

(3) the eigenvalues form the geometric series 1, 1/a, 1/a 2 , . . . , l/a n ~ l ; 

(4) the right eigen vectors are independent of a and have the simple form: 
V m (i) = (i- l)^ 1 ^ 1 ) + (-l) n - l+m (T-) f° r V« m , m > 1. 

Proof. The cross-symmetry (1) follows from Proposition l3.14 and the multiplicative property 
(2) follows from the shuffling interpretation and equation ([2]). Property (1) implies that the 
eigen structure is quite constrained. Properties (3) and (4) follow from results of Cuicu 

a- □ 

The following Corollary also follows as a special case of Theorem 2.2 in [5]. 

Corollary 3.3. Consider a deck of n cards with the ace of spades starting at the bottom. 
The chance that the ace of spades is at position j from the top after an a-shuffte is 

(6) Q a (j) =Pa(n,j) = ±J2(k-ir~jki- 1 . 

a ti 



From the explicit formula, we are able to give exact numerical calculations and sharp 
asymptotics for any of the distances to uniformity. The results below show that log 2 n + c 
shuffles are necessary and sufficient for both separation and total variation (and there is a 
cutoff for these). This is surprising since, on the full permutation group, separation requires 
2 log 2 n + c steps whereas total variation requires | log 2 n + c. Of course, for any specific n, 
these asymptotic results are just indicative. 

Remarks on Table [0 We use Proposition 13.11 to give exact results when n = 52. For 
comparison, the upper table gives exact results for the full deck using [3]. The lower table 
shows that it takes about half as many shuffles to achieve a given degree of mixing for a 
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Table 1, Distance to uniformity for a deck of 52 cards. The upper table 
assumes distinct cards, and the lower table follows a single card starting at 
the bottom of the deck. 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


TV 


1.00 


1.00 


1.00 


1.00 


.924 


.614 


.334 


.167 


.085 


.043 


.021 


.010 


SEP 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


.996 


.931 


.732 


.479 


.278 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


TV 


.873 


.752 


.577 


.367 


.200 


.103 


.052 


.026 


.013 


.007 


.003 


.002 


SEP 


1.00 


1.00 


.993 


.875 


.605 


.353 


.190 


.098 


.050 


.025 


.013 


.006 



card at the bottom of the deck. For example, the widely cited '7 shuffles' for total variation 
drops this distance to .334 for the full ordering, but this requires only 4 shuffles to achieve 
a similar degree of randomness for a single card at the bottom. 



For asymptotic results, we first derive an approximation to separation, which also serves 
as an upper bound for total variation. Finally, we derive a matching lower bound for total 
variation. Proofs have been omitted for brevity, but again full details are available in [2]. 

Proposition 3.4. After an a-shuffle, the probability that the bottom card is at position i 
satisfies 

1 a n ~ i+l „ , 1 



< Qaii) < - 



a 1 — a n ~"" w al — a n 1 ' 
where for brevity we have set a = 1 — 1/a. In particular, the separation distance satisfies 

1 _ It" < SEP(a) < 1 - -- 

a 1 — a n a 1 



If a = 2 log2 ( n ) +c = n2 c , then our result shows that the SEP(a) is approximately 

1 



1 e- 2 ~ 



2 C 1 - e - 2_c ' 

and for large c this is ~ 2 -c ~ 1 . The fit to the data in Table Q] is excellent: for example after 

26 
1021 



ten shuffles of a fifty-two card deck we have 2 c 1 = which is very nearly the observed 



separation distance of 0.025. 

Remark 3.5. Proposition \3.4\ gives a local limit for the probability that the original bottom 
card is at position j from the bottom. When the number of shuffles is log 2 n+c, the density of 
this (with respect to the uniform measure) is asymptotically z(c)e~^ 2C , with z a normalizing 
constant (z{c) = l/2 c (e J//2C — 1) ). The result is uniform in j for c fixed, n large. 

Proposition 3.6. Consider a deck of n cards with the ace of spades at the bottom. With 
a = 1 — 1/a, the total variation distance for the mixing of the ace of spades after an a-shuffle 
is at most 

a n+1 aa 2 (l-a n - v ) 1 , f al 

H log 

1 — a n n(l — a n ) ralog(l/a) \n c 
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and at least 

a n a{l-a n ) 1 fal-a^ 1 



H log — 

a ra_1 na(l — a n_1 ) nlog(l/a) Vn a"" 1 



After log 2 n + c shuffles, that is when a = 2 c n, Proposition 13,61 shows that the total 
variation distance is approximately (with C = 2 C ) 

i/c ,,\ . l-Clog(e»/c_l) 



Clog^-l)j + (el/c _ i} . 

Thus when c is 'large and negative,' the total variation is close to 1, and when c is large 
and positive, the total variation is close to 0. Thus total variation and separation converge 
at the same rate. This is an asymptotic result and, for example, Table Q] supports this. 

Similar, but more demanding, calculations show that if the ace of spades starts at position 
i, and max(i/n, (n — i)/n) > A > for some fixed positive A, then ^\og 2 n shuffles suffice 
for convergence in any of the metrics. We omit further details. 

4. Separation distance for the general case 

A main result of Bayer and Diaconis [3] is the simple formula for an a-shuffle of a deck 
of n distinct cards: 

(7) aw4(" + * 

where r = r(cr) is the number of rising sequences in a, equivalently one more than the 
number of descents in a . This formula allows simple closed form expressions for a variety 
of distances as well as asymptotic analysis. 

In this section we work with general decks containing Di cards labelled i, 1 < i < m. 
The formulae of this section hardly resemble the elegant expression above. Further, we only 
give precise formula for the least likely deck. The following lemma shows that this deck, 
where the separation distance is achieved, is the reverse the initial deck configuration. This 
is equivalent to Theorem 2.1 from [5]. 

Proposition 4.1. Let D be a deck as above. After an a-shuffle of the deck with 1 's on top 
down to m 's on bottom, the least likely configuration is the reverse deck w* with m 's on top 
down to 1 's on the bottom. 

Proof. The only cuts of the initial deck resulting in w* are those containing no pile with 
distinct letters. For all such cuts, each rearrangement of the deck is equally likely to 
occur. □ 

While finding a completely general formula for Q a {w) for arbitrary w is infeasible, below 
we do this for w* . 

Theorem 4.2. Consider a deck with n cards and Di cards labeled i, i = 1, . . . , m. Then the 
separation distance after an a-shuffle of the sorted deck (1 's followed by 2 's, etc) is given by 

SEP(a) = l--f U D \J2 (a-^O^nC^-Vi^-fe-Vi-l)^)- 

1 • • • mj - =l 

0=ko<--<k m -i<a 
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Proof. Prom the analysis in the proof of Proposition 14. 1\ Q a (w*) is given by 

Qa(w*)= J2 -^( A A )- ( n V' 

Al+ ... +Aa=n a \Al,...,A a J { Dlr __ iD J 

A refines D 

where l A refines means there exist indices ki, . . . , k m _\ such that A\ + • • • + = D\ 
and, for % = 2, ... ,m — 1, A^^+i + • • • + A^ = D{. Taking the fcj's to be minimal, the 
expression for Q a (w*) simplifies to 

^ m— 1 

(8) - 52 («-^-o Dm n ((%-%-i) Dj - • 

0=fco<---<fe m _i<a j=l 

The result now follows from Proposition 14,11 □ 



Table 2. Separation distance for k shuffles of 52 cards. 



k 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


BD-92 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


.995 


.928 


.729 


.478 


.278 


blackjack 


1.00 


1.00 


1.00 


1.00 


.999 


.970 
















1.00 


.997 


.997 


.976 


.884 


.683 


.447 


.260 


.140 


.073 








1.00 


1.00 


.993 


.875 


.605 


.353 


.190 


.098 


.050 


.025 


.013 


.006 


redblack 


.890 


.890 


.849 


.708 


.508 


.317 


.179 


.095 


.049 


.025 


.013 


.006 


mm 


1.00 


1.00 


.993 


.943 


.778 


.536 


.321 


.177 











Remarks on TablelM We calculate SEP after repeated 2-shuffles for various decks using The- 
orem H2J (blackjack) 9 ranks with 4 cards each and another rank with 16 cards; (JfcO'v'^) 4 
distinct suits of 13 cards each; (A<|b) the ace of spades and 51 other cards; (redblack) a two 
color deck with 26 of either color; and ([o|+[¥[5]j^|) a deck with 5 cards in each of 5 suits. 
The missing entries in Table [2] highlight the limitations of exact calculations using Theorem 
PI 



Remark 4.3. Comparing the data in Tabled for A<|b and redblack shows that these two 
cases are remarkably similar. Indeed, both cases exhibit the same asymptotic behavior, which 
is remarkable since the A<(k has a state space of size 52 while redblack has a state space of 
size around 5 x 10 14 . 

Now we derive a basic asymptotic tool which allows asymptotic approximations for gen- 
eral decks. 

Proposition 4.4. Let m > 2 and a be natural numbers, let £1, . . ., £ m be real numbers in 
[0,1]. Let 

r i; ■ ■ ■> r m be natural numbers all at least r > 2. Let 

S m (a;£,r)= J] (a x + ^) ri • • • (a m + UT m - 

ai,...,a m >0 
a 1 +...+a m =a 



s 
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Then 



S m (a; £,r) - n!---^! +& + ... + ^i+.-.+^+m-t 

(ri + . . . + r m + m - 1)1 

.V^V^-iV 1 y>+£i + ---+£m) ri+ --- +rm+m - 1 - 2j 



j J \3(r - I) J (n + . . . + r m + m - 1 - 2j)! 



Consider a general deck of n cards with Di cards labelled i. We use Proposition 14.41 to 
find asymptotics for the separation distance given in Theorem 14.21 The following is our 
'rule of thumb.' 

Theorem 4.5. For a deck of n cards as above, suppose Di > d > 3 for all 1 < i < m. 
Then we have 

m — 1 m— 1 / -. \ -ii 

a r—k. ;(m-l\( j\n+m-l 



SEP(a) = 1 - (1 + r,)- f — — f Y ( m . 

KJ u (n + l)---(n + m-l) ; \ j J 

where n is a real number satisfying 

/ n 2 \m-l 

N<(1+ ^ 7^Z II I Tv? ) "I 



1 - - 



3(d-2)(a-m + l) 2 . 
Proof. To evaluate the expression in Theorem 14.21 we require an understanding of 

m— 1 „\ p\ m— 1 

E ^ ll(«f J -(oi-i) J, o = / ■••/ E II (^(ai-i-^)"'- 

oi+.. .+a m =a j=l ai+...+a m =a jf=l 

a j > 1 aj > 1 

We now invoke Proposition 14.41 Thus the above equals for some \6\ < 1 
j^.i Z 1 ... f 1 / (a-Cm-lj + gx + .-. + g^-!)" | 



m— 1 



o Jo v ni 
m-lV 1 y(a-(m-l)+£i + ...+£m-iV^ 2/ 



We may simplify the above as 

and evaluating the integrals above this is 

,2 \ m-l \ n , _ _ _ n , m-l 



a-j) n - m+1 . 



j=0 

The Theorem follows. □ 

For simplicity we have restricted ourselves to the case when each pile has at least three 
cards. With more effort we could extend the analysis to include doubleton piles. The case 
of some singleton piles needs some modifications to our formula, but this variant can also 
be worked out. Below we use our rule of thumb to calculate separation for the same decks 
as in Table [2j 
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Table 3. Rule of Thumb for the separation distance for k shuffles of 52 cards. 



k 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


BD-92 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


.995 


.928 


.729 


.478 


.278 


blackjack 


1.00 


1.00 


1.00 


1.00 


.999 


.970 


.834 


.596 


.366 


.204 


.108 


.056 




1.00 


1.00 


.997 


.976 


.884 


.683 


.447 


.260 


.140 


.073 


.037 


.019 


redblack 


.962 


.925 


.849 


.708 


.508 


.317 


.179 


.095 


.049 


.025 


.013 


.006 




1.00 


1.00 


.993 


.943 


.778 


.536 


.321 


.177 


.093 


.048 


.024 


.012 



Remarks on Tabled The first row gives exact results from the Bayer-Diaconis formula for 
the full permutation group. The other numbers are from the rule of thumb. Roughly, the 
single card or red-black numbers suggest that half the usual number of shuffles suffice. The 
Black-Jack (equivalently Baccarat) numbers suggest a savings of two or three shuffles, and 
the suit numbers lie in between. The final row is the rule of thumb for the Zener deck with 
25 cards, 5 cards for each of 5 suits. 

While asymptotic, Theorem 14.51 is astonishingly accurate for decks of practical interest. 
For instance, comparing exact calculations in Table[2]with approximations using this rule of 
thumb in Table [3] shows that after only 3 shuffles, the numbers agree to the given precision. 
Moreover, the simplicity of the formula in Theorem 14.51 allows much further computations 
than are possible using the formula in Theorem 14.21 

We now give a heuristic for why our rule of thumb is numerically so accurate. For k > 0, 
define 



fk{z) = jr 



r k z r _ Mz) 



where Ak{z) denotes the A;-th Eulerian polynomial. The sum over a±, . . ., a m appearing 
in our proof of Theorem 14.51 is simply the coefficient of z a in the generating function (1 — 
z yn-i f Di (z) ■ ■ ■ fn m (z). Our rule of thumb may be interpreted as saying that 

(9) (1 - zT- l f Dl (z) ■ ■ ■ f Dm (z) » P 1 ^" ^ , (1 - zT^f n+m ^(z). 

[n + m— lj! 

To explain the sense in which ([9]) holds, note that fk(z) extends meromorphically to the 
complex plane, and it has a pole of order k + 1 at z = 1. Moreover it is easy to see that 
fk(z) — k\/(l — z) k+1 has a pole of order at most k at z = 1. Therefore, the LHS and 
RHS of ([9j) have poles of order n + 1 at z = 1, and their leading order contributions match. 
Therefore the difference between the RHS and LHS of ([9]) has a pole of order at most n at 
z = 1. But in fact, this difference can have a pole of order at most n — d at z = 1, and 
thus the approximation in ([9]) is tighter than what may be expected a priori. To obtain our 
result on the order of the pole, we record that one can show 

v.] / ( z — I) \ k+l 

'»M = (r^BT(W) +«-*>+ °a 
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5. GlLBREATH PRINCIPLE AT WORK 

Conger and Viswanath note that the initial configuration can affect the speed of conver- 
gence to stationary. Perhaps this is most striking in the case of Section [3] where a single 
card is tracked. Recall Table [H giving calculations for the distinguished card beginning at 
the bottom of a deck of 52 cards. In contrast, Table H] gives calculations for the distin- 
guished card starting in the middle, at position 26. For the latter, both total variation and 
separation are indistinguishable from zero after only four shuffles. 



Table 4. Distance to uniformity for a single card starting at the middle of 
a 52 card deck. 





1 


2 


3 


4 


TV 


.494 


.152 


.001 


.000 


SEP 


1.00 


.487 


.003 


.000 



Consider next a deck with n red and n black cards. First take the starting condition of 
all reds atop all blacks. If the initial cut is at n (the most likely value) then the red-black 
pattern is perfectly mixed after a single shuffle. More generally, the chance of the deck w 
resulting from a single 2-shuffle of a deck with n red cards atop n black cards is given by 



Q2(«0 = 2sr(2 h "" ) + 2 11 "-' 



i 



where h(w) is the number of red cards before the first black card and t(u>) is the number of 
black cards after the final red card; see [2]. In particular, the total variation after a single 
2-shuffle is 
(10) 
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Evaluating this formula for 2n = 52 give a total variation of 0.579. 

Now take the starting condition to alternate red black red black, etc. As motivation, we 
recall a popular card trick: Begin with a deck of 2n cards arranged alternately red, black, 
red, black, etc. The deck may be cut any number of times. Have the deck turned face up 
and cut (with cuts completed) until one of the cuts results in the two piles having cards of 
opposite color uppermost. At this point, ask one of the participants to riffle shuffle the two 
piles together. The resulting arrangement has the top two cards containing one red and one 
black, the next two cards containing one red and one black, and so on throughout the deck. 
This trick is called the Gilbreath Principle after its inventor, the mathematician Norman 
Gilbreath. It is developed, with many variations, in Chapter 4 of [10]. From the trick we 
see that beginning with an alternating deck severely limits the possibilities. Analyzing the 
trick reveals the following formula, 



(11) 



-,2n 



Q 2 H 



2 n- 
2 r 



if w is the initial alternating deck, 
if w can result from an odd cut, 
if w can result from an even cut, 
otherwise, 



RIFFLE SHUFFLES WITH REPEATED CARDS 
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where an odd (resp. even) cut refers to the parity of cards in either pile. From this we 
compute 

1 / 2 n 4- 2 n_1 - 1 \ 
(12) ll32-f/|| TV = Ml-^^ J, 

which goes to .5 exponentially fast as n goes to infinity, and indeed is already .500 for 
2n = 52. In contrast, starting with reds above blacks, asymptotic analysis of (|10p shows 
that the total variation tends to 1 after a single shuffle when n is large. Thus again an 
alternating start leads to faster mixing. 
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